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O I Abstract. New results on transverse mass spectra of neutral pions measured at central rapidity are presented for impact 

^ . parameter selected 200A-GeV S + S and S + Au collisions. The distributions cover more than 8 orders of magnitude 

C ' in cross section over the range 0.3 GeV/c < Pt < 4.0 GeV/c. Detailed comparisons to results from pp collisions 

^ I are made. The spectra from all systems show a clear power-law like shape with similar curvature. Collisions of S + Au 

• 1^ . exhibit a larger mean transverse momentum than pp increasing with centrality. Predictions of string models and by 

' hydrodynamic approaches including collective expansion and decays of short lived resonances are compared to the 

' data and the implications are discussed. 



1 Introduction 

Heavy ion reactions at sufficiently high energies are able to 
create a strongly compressed and highly excited reaction zone 
which may provide the conditions for the transition to a decon- 
fined state of nuclear matter, the Quark-Gluon Plasma. Direct 
information on this state and the initial conditions are provided 
by probes which decouple early in the reaction, like photons 
and dileptons. Hadrons which interact strongly decouple much 
later and their transverse momentum spectrum therefore gives 
information about the condition of the system at freeze out. By 
combining the information from photons or dileptons and from 
hadrons and comparing with data from proton-proton interac- 
tions it is hoped that the development of the heavy-ion system 
from the initial state of high density and temperature through 
freeze out can be reconstructed. In early pp scattering experi- 
ments, it has been found that px spectra of pions in the range 
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up to about 2 GeV/c exhibit a nearly exponential shape with a 
slope constant of about 160 MeV/c. This is at first sight sur- 
prising since most of the pions are not produced directly but 
emerge from heavier resonances with different decay energies. 
Furthermore, exponential spectra are usually observed in sys- 
tems, where thermodynamical equilibrium is reached through 
a sizable number of rescatterings. In a proton-proton collision 
this is hardly conceivable. Using the bootstrap hypothesis and 
particular assumptions on the particle mass spectrum, Hage- 
dorn was able to describe the exponential shape of the spectra 
and the limiting value of the slopes [|l|] . 

Proton-proton experiments at the CERN ISR extended mea- 
surements to Pt values beyond 2 GeV/c [^[^|| and subse- 
quently revealed distinct deviations from an exponential shape 
approaching a more powerlaw-like behavior This was inter- 
preted as the onset of "hard" parton scattering [||] and can be 
treated to some extent by perturbative QCD [^]. 

Data on hadron-nucleus collisions display further charac- 
teristic differences when compared to proton-proton collisions. 
These are often parametrized by a factor A°'^'P'^^ connecting 
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p+p and p+A cross sections. At low transverse momenta, a{pT) 
turns out to be about 0.7, which is close to the value of 2/3 ex- 
pected from a geometrical picture of hadron production from 
an opaque target nucleus. For larger transverse momenta, a{pT) 
increases, which is interpreted as the onset of parton scatter- 
ing extending over the whole nuclear volume. In this picture, 
a{pT) should finally approach the value 1. The surprising ob- 
servation of a > 1 in the experimental spectra at high pt ^ ^ 
GeV/c (Cronin-effect [@]) was then interpreted as the effect of 
multiple parton scattering [|^,^. 

All of the features of proton-proton and hadron-nucleus 
scattering described above are expected to be present in en- 
ergetic collisions between heavy nuclei. In addition, due to the 
nuclear matter being highly excited and compressed over a large 
volume, new specific thermodynamic or hydrodynamic features 
are expected to become apparent. 

Theoretically, both the description of nucleus-nucleus colli- 
sions as a superposition of hadron-hadron collisions [[l^,|ll]Jl^ 
as well as the assumption of a hot hydrodynamically expand- 
ing system of hadrons or partons [ pjHl4| ], originally proposed 
by Landau JI^ are being pursued. 

In both approaches, the approximate exponential shape of 
the final spectra at low px results from a convolution or super- 
position of various basic processes and sequential decays, and 
cannot easily be associated with a "temperature". At high pT, 
on the other hand, it is expected that the spectra converge to the 
"hard scattering" shapes observed in proton-proton scattering 
and successfully described by perturbative QCD. 

In order to distinguish between the theoretical approaches 
and to isolate possible new features of heavy-ion collisions, it 
is mandatory to investigate the transverse momentum spectra 
of identified particles over a large range in px extending in par- 
ticular to the hard scattering regime. 

In general, the invariant differential cross sections of pro- 
duced particles measured from a given reaction system depend 
not only on the transverse momentum pT, but also on the avail- 
able center-of-mass (cms) energy -/i, the rapidity y of the par- 
ticle, and the associated centrality. A comparison of pr spectra 
from different collision systems or different data sets therefore 
has to carefully consider the dependencies on the quantities 
noted above. 

In JI^, preliminary results of neutral pion spectra for S- 
induced reactions have been presented. In the present paper, 
the final results on high precision neutral pion px spectra from 
collisions of 200 ^-GeV sulphur nuclei (y/s = 19.4 AGeV) 
with sulphur and gold targets in the range 0.3 GeV/c < pr < 
4.0 GeV/c and 2.1 < y < 2.9 for various centralities after a 
complete reanalysis are shown. They will be compared to data 
from p H- p and p H- A interactions. The production of 77 mesons 
measured within the same experiment over the px range of 
0.5 < ^ 3.5 GeV/c [ |l7t and an upper limit on direct photon 
production i jlsj l have been recently published. 

2 Experiment 

The WA80 experiment at the CERN-SPS (Fig. |l]) was equipped 
with a finely segmented electromagnetic spectrometer com- 
posed of 3798 lead-glass modules with photomultiplier read- 
out. The lead-glass was arranged into three independently cal- 
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Fig. 1. WA80 experimental setup 



ibrated arrays, of roughly similar size. Two arrays consisted of 
TFl lead-glass of 4.0 x 4.0 x 40 cm^ (15 Xq) deployed as 
towers to the left and right of the beam axis. The third array, lo- 
cated below the beam axis, was the S APHIR lead-glass detector 
[ pO| ] already used in the WA80 ^^O run period [ pll and which 
consisted of SF5 lead-glass modules of 3.5 x 3.5 x 46 cm^ 
(18 Xq). The entire photon spectrometer, which was located 
at a distance of 9 m downstream of the target, provided cov- 
erage of 10% to 50% of full azimuth 4> over the polar range 
6.2° < i3iab < 13.9°. This resulted in a rather uniform accep- 
tance in the rapidity range 2.1 < y < 2.9 for reconstruction of 
7r° and 77 mesons via their 77 decay mode. 

Immediately in front of the photon spectrometer was a 
double-layer charged particle veto (CPV) counter which cov- 
ered the lead-glass region of acceptance. Each layer of the CPV 
consisted of streamer tubes with charge-sensitive pad readout, 
with pads of dimension similar to the lead-glass modules [|2|. 

The minimum bias (mb) trigger requires a valid signal of 
the beam counters and a minimum amount of transverse en- 
ergy Et ^ 5GeV, detected by the Mid-Rapidity Calorime- 
ter MIRAC [HI. Data have been taken with the 200 AGeV 
sulfur beam on targets of Au (250 mg/cm^) and S (205 and 
510 mg/cm^). For the present analysis 7.9 million Sh-Au and 
1.4 million Sh-S minimum bias events were accumulated, al- 
lowing precise reconstruction of the tt" pr distributions over a 
wide range of 0.3 < pr < 4.3 GeV/c for S + Au reactions. The 
minimum bias cross section has been calculated from the num- 
ber of beam and minimum bias triggers and the target thick- 
nesses, and was corrected for small contributions of background 
effects (typically a few percent). As a result, we obtain (Tnif) = 
1450 mb and 3600 mb for Sh-S and Sh-Au reactions, respec- 
tively. The S + Au value is approximately 20 % larger than for 
the 1987 WA80 data Q, due to the more biased trigger set- 
ting used then. These absolute cross sections have an overall 
systematic error of 10%. 
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3 Pion Reconstruction and Efficiency 

For the extraction of the neutral pion yield via the 77 decay 
branch, all showers in the lead-glass have been considered. 
While additional criteria like a lateral dispersion cut [ [l9[ ] or 
the use of the CPV can improve the quality of photon identi- 
fication, these cuts are of little importance for the analysis of 
neutral pions which are uniquely identified via their invariant 
mass peaks. 

Hits in the detector are combined to pairs to provide dis- 
tributions of pair mass vs. pair transverse momentum for all 
possible combinations. These distribution are obtained both 
for real events (R{minv,PT)) and for so-called mixed events 
(M{minv ,Pt)), where a hit from one event is combined with a 
hit from another event with similar multiplicity. M{minv,PT) 
provides a good description of the combinatorial background. 
It is subtracted from R{minv,PT) to obtain the mass distribu- 
tion of neutral pions. 



class 


Et (GeV) 


<J /o"mb 


cell 
occupancy 


no. of 
events 


1 


< 20.0 


30.0 % 


0.17 % 


2.37 • 10" 


2 


20.0- 42.0 


20.1 % 


0.43 % 


942 • 10^ 


3 


42.0- 66.3 


15.2 % 


0.82 % 


461 • 10^ 


4 


66.3- 86.0 


9.5 % 


1.24% 


539 • 10^ 


5 


86.0-107.5 


10.5 % 


1.55 % 


1.48 • 10" 


6 


107.5-121.5 


7.0% 


1.90% 


1.02 • 10" 


7 


121.5-133.5 


4.9 % 


2.10% 


717- 10^ 


8 


>133.5 


2.8 % 


2.40 % 


407 • 10^ 



Table 1. Centrality classes, as selected by the amount of transverse 
energy observed in MIRAC, and cell occupancies for ^^S + Au colli- 
sions. 



In central reactions of ^^S + Au the average particle oc- 
cupancy in the detector reaches values of 2.4 %/module (see 
Table p. This leads to a finite probability that showers in the 
detector overlap and influence each other. Such an overlap may 
result in the following modifications of the measured signals: 

1 . The signal of one shower is "absorbed" by another one - 
this happens mostly to low energy showers. In such a case, 
the photon is lost and if it originated from a vr" -decay, the 
tt" parent particle itself cannot be reconstructed and is lost. 

2. The signal of one shower is influenced by that of a nearby 
shower. In this case, the measured energy and/or position 
may be changed and, thereby, the momentum of the pho- 
ton may be changed. If it originated from a neutral pion, 
depending on how important these changes are, this may 
either result in 

(a) changes of the apparent invariant mass such that the 
pion cannot be reconstructed, or 

(b) changes that still allow to identify a pion by its proper 
invariant mass, but at a different momentum. 

In addition, case |b| leads to a modification of the shape of 
the invariant mass distribution of 7r"s, which has to be taken 
into account in the extraction of the peak content. 



These effects lead to a detector efficiency for pion recon- 
struction which depends on particle density. To study this de- 
tection efficiency, the GEANT [|5]| simulation package has been 
used to create artificial signals for the lead glass modules corre- 
sponding to neutral pions hitting the detector In total, w 10^ 7r°s 
have been simulated with uniformly distributed transverse mo- 
menta. 
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Fig. 2. Geometrical acceptance for neutral pions as a function of pr- 



The signals were then analyzed with the same analysis pro- 
grams as used for real data. Analyzing them as single particles, 
allows to study detector effects and the analysis programs. To 
obtain a tt" efficiency, however, these simulated photon pairs 
have to be studied in a realistic particle density environment. 
This has been done by superimposing them onto real measured 
events. Effects of detector noise and digitization of the photo- 
multiplier signals are implemented. Most importantly, the sim- 
ulation provides the means to extract the probability for a pion 
at given input transverse momentum p^-* to be measured with 

p^\ These probabilities measured as a function of can be 
understood as a matrix which transforms real physical distribu- 
tions into measured ones. The efficiency matrix contains both 
effects of detection efficiency as described above, as well as 
the geometrical acceptance. The acceptance can be easily sep- 
arated and is shown in Fig.^ as a function of the transverse 
momentum. The efficiency matrix is calculated for each event 
class (i.e. given particle density). The representation as a ma- 
trix has the advantage that the actual shape of the momentum 
distributions does not need to be known beforehand. 

The efficiency matrices are used to correct the measured 
distributions in the following way: 



A reasonable tt'^ distribution is chosen as initial ansatz. 
The distribution is transformed by applying the efficiency 
matrix. 

The resulting distribution (B) is compared to the actual 
measured distribution (A). 

The input distribution is modified by multiplying it with the 
ratio A/B. 

Steps 1-4 are iterated until A and B agree. 
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The result of this procedure is shown in Fig. ||, where the 7r° 
reconstruction efficiency[] is plotted for the most peripheral and 
most central S + Au reactions, respectively. It can be seen that 
for very peripheral reactions, the correction is relatively small, 
particularly for > 1 GeV/c. However, for the most central 
reactions, i.e. for the highest particle densities, the correction 
becomes much more important. The general behavior can be 
understood from the previous discussion. The drop in the 7r° 
reconstruction efficiency at low pr is largely caused by the in- 
creasing probability with decreasing energy to misidentify or 
lose a low energy photon in the radial tail of a higher energy 
shower, such that the (mostly low pt) neutral pion cannot be 
reconstructed and is lost for further analysis. At large trans- 
verse momenta, the reconstruction efficiency exceeds the value 
1 . This is an effect of limited transverse momentum resolution 
on steeply falling spectra; effectively the observed distribution 
is broadened over the original one by the apparatus, thus giving 
rise to 6-^0 > 1 (item 2b in the list above). 

This is illustrated in Fig.Q Here, pions have been simu- 
lated with transverse momenta in a narrow range (0.8 GeV/c < 
Pt < 1-0 GeV/c) The resulting pT distributions after the over- 
lap procedure discussed above are shown as dark grey his- 
tograms. One can see that already in peripheral collisions (up- 
per right) there is some broadening of the distribution, the ef- 
fect is however much stronger in central collisions (lower right). 
Here a significant fraction of the pions will be measured with 
higher transverse momentum. For the steeply falling spectra 
indicated as line histograms this may result in a larger appar- 
ent number of pions in the higher pT bins and thus lead to an 
"efficiency" > 1. An equivalent broadening effect can be seen 
in the experimental data e.g. in the invariant mass spectra dis- 
played on the left hand side of Fig. ^ Overlap effects lead to an 
additional high mass tail of the pion peak in central collisions. 

The systematic errors on the measured transverse momen- 
tum spectra are dominated by the following contributions^: 



We refer now to the "pure" efficiency without any effect from the 
geometrical acceptance. 

^ The estimates given here apply to central collisions where the er- 
rors are largest. 
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Fig. 3. Neutral pion reconstruction efficiency for very peripheral (a) 
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Fig. 4. Influence of shower overlap on the neutral pions with 
0.8GeV/c < Pt < l.OGeV/c for peripheral and central colli- 
sions. The plots on the left hand side show two-photon invariant mass 
spectra, where the light grey area indicates the integration range for 
the pion extraction. The dark grey histograms on the right hand side 
show the Pt distribution of pions simulated with 0.8 GeV/c <Pt< 
1.0 GeV/c after overlap effects from other showers are included. For 
reference, the relevant region of the pT-spectra is also shown. See text 
for further discussion. 



- An uncertainty in the absolute calibration of the momentum 
scale of 1%. This may be translated into an uncertainty in 
the yields| of a few % at low px rising to « 13% at px = 
3.5 GeV/c. 

- Uncertainties in the 7r° extraction which include the error 
in the determination of the invariant mass peak content and 
the propagation of the energy resolution through the accep- 
tance and efficiency corrections. This leads to an error of 
16% at Pt — 0.5 GeV/c which decreases to below 3% at 
highpT- 

These two major contributions to the systematic errors and their 
quadratic sum are shown in Fig. ||. 
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Fig. 5. Systematic errors of the extraction of the tt" yields for central 
S + Au collisions as a function of pT. 



and very central S + Au reactions (b). 



^ This error contribution is much smaller in the ratios of photon or 
pion spectra as e.g. used for the direct photon analysis mM- 
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4 Results 

4.1 Nuclear Reactions 

Efficiency corrected tt" invariant cross sections measured in 
the rapidity range 2.1 < y < 2.9 are presented in Fig.|^ for 

5 + S and S + Au reactions. The data are shown for central, 
peripheral and minimum bias trigger conditions. The central 
event class of 7.7% of (Tmb in S + Au collisions corresponds 
to complete overlap of the S nucleus with the Au target, with 
an average of w 110 participating nucleons, to be compared to 
an average of sa 11 participating nucleons for the 30% most 
peripheral events. From a first inspection it is seen that the 
slopes appear systematically flatter for central collisions than 
for peripheral ones. This might indicate a higher temperature 
in more violent central collisions. Although the slope parame- 
ters cannot be identified a priori with a temperature in a ther- 
modynamic sense, we shall often use the label "temperature" 
in the following.^ Superimposed on the data are exponentials 
da/dm\, oc exp(— tot/T) fitted to the region 0.8 GeV/c^ 
< rriT < 2.0 GeV/c^. The fitted slope parameters are indi- 
cated. In this parameterization of the intermediate my data, the 
flattening of the distributions with increasing centrality results 
in a higher apparent temperature. The high precision and the 
large rriT range of the present 7r° data allows to recognize, how- 
ever, that the spectral shapes are only very poorly described by 
an exponential curve. Instead, a concave behavior reminiscent 
of pp data at similar or higher y^s is observed over the whole 
rriT range. Its properties are more systematically analyzed by 
calculating the local slope, defined as 



diriT 



log ( 



dp 



(1) 



The results are plotted in Fig. ^ for the S + Au data. Here, each 
individual value of Ti^c has been determined locally from three 
adjacent data points of Fig. ^ Such an analysis puts very high 
demands on the statistical and systematic uncertainties in or- 
der to keep point-to-point fluctuations at an acceptable level. 
Shown together with the heavy ion data are pp minimum bias 
data which will be discussed in section 4.2. 

As expected, the observed local slopes are not constant, 
as would be the case for a purely exponential distribution, but 
change continuously over the measured rriT range both for pe- 
ripheral and central collisions. This again indicates that fits 
with one or two exponentials, which have often been presented 
elsewhere (e.g. Ref. [|6|]), can only provide a crude approxima- 
tion valid over a restricted range. Fig.^ also demonstrates 
that the slope parameters in the intermediate rriT range are sys- 
tematically higher for central compared to peripheral data. 

In order to systematically compare data from different sys- 
tems and/or taken under different trigger conditions it is often 
more convenient to apply a description in terms of a functional 
form. According to the above discussion, such a parameteriza- 
tion must allow for pr dependent variations in the curvature of 
the spectra. At high pr, where perturbative QCD becomes ap- 
plicable, the spectra are expected to attain a power-law behav- 
ior, as observed in numerous high energy pp measurements (see 

All slope parameters will be given in units MeV. 



for example Ref. [pT]]). The heavy-ion data of this experiment 
seem to follow that trend even into the intermediate px range. 
Therefore, a parameterization originally inspired by QCD 
and successfully applied already to pp data has been cho- 
sen to fit to the spectra: 



dp^ 



Po 



Pt +Po 



(2) 



with C, Pa and n taken as free parameters. A link to the more 
familiar exponential slope parameter T is obtained from the 
derivative of Eq.0 according to 



power — law 



Pa 



1 



■PT 



(3) 



Thus, po/n characterizes the slope of the transverse mo- 
mentum spectra in the limit px 0, while 1/n characterizes 
its gradient along px, i.e. the strength of the concave curvature. 
The analytical form of Eq.|^ can in that way be compared di- 
rectly to the numerical values of the local slopes presented in 
Fig.0. As stressed by Hagedorn [p8|], and immediately recog- 
nized from Eq. ^ po and n are highly correlated. We have there- 
fore chosen to express Eq. ^ as a function of the local slopes Ti 
and T2 at two fixed values of px (IGeV/c and 3GeV/c) - 
these are less correlated and therefore better suited as fit pa- 
rameters. 

Fits to the data using Eq. || are presented in Fig. ^. The ex- 
tracted parameters are compiled in Table ^ together with the 
results of analogous fits to pH-p data [|^]. All slope parame- 
ters are observed to increase smoothly with the centrality of 
the reaction, e.g. Ti increases from approximately 196 MeV in 
peripheral collisions to approximately 206 MeV in central col- 
lisions. 

Another way to investigate the differences between periph- 
eral and central data in a model independent way is to study the 
ratio of the transverse momentum distributions as a function of 
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Fig. 7. Local slope presentation of central and peripheral S+Au tt" 
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data drawn together with minimum bias pp data 
line shows the local slopes calculated according to Eq. 3| from the fit 
of Eq.^to the pp data. 
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Fig. 6. Invariant cross sections of tt" mesons from reactions of S + S (left) and S + Au (right) at 200j4-GeV measured in the rapidity range 
2-1 < 2/ < 2.9. The events are selected for centrality with the percentage of the minimum bias cross section as indicated. The lines drawn with 
the data are exponentials fitted to the region 0.8 GeV/c^ < rriT < 2.0 GeV/c^ with slope parameters as indicated. 
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(30%) S + Au reactions as a function of pr- 



Pt- Fig. ^ shows that with increasing pr, the pion multiplicity 
in central data is enhanced relative to that of peripheral data. 
Measured over the dynamic range of the S + Au data, a varia- 
tion by a factor of > 2 is observed. Similar findings in the cross 



section ratios have been reported earlier by us for central and 
peripheral O + Au reactions at 200 GeV [Q. 

The mean transverse momentum (pt) is an observable con- 
taining important information about the reaction dynamics of 
the system and is often used to provide a convenient means of 
comparison among different experiments and model calcula- 
tions. The evaluation of (px) is usually achieved analytically 
by using parameterizations of the experimental distributions 
(see below). Unless the cross section is measured down to very 
low Pt and the functional form describes the data well, this 
method may lead to systematic errors due to the dominance of 
the low Pt cross section. If only the variation of (pt) as a func- 
tion of another observable, like for example the centrality, is of 
importance, it is more advantageous to calculate instead the 
truncated mean transverse momentum above a certain thresh- 
old directly from the experimental data. This can be done by 
the relation: 



'T d 



' dN , 

PT-: dpT 

dpT I 



dN 
dpT 



dpi 



(4) 



Here, d — 0.4 GeV/c is chosen according to the lowest pt 
data point where systematic uncertainties imposed by accep- 
tance and efficiency corrections are still well under control. 
For a purely exponential invariant cross section, i.e. da/dp^ cx 
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Fig. 8. Fits to 200 A GeV S + S and S + Au data using Eq.^ with parameters given in Table ^. 



exp{—pT/T), and for d — 0, one immediately finds {pT)d=o — 
2T. For finite values of d, however, this simple relation does not 
hold and {pT)d decreases as a function of d. 

Most interesting in connection with the thermodynamical 
picture of the compressed and excited heavy ion system is the 
variation of (px) with Et or the multiplicity dNch/di] of pro- 
duced particles. We therefore have studied {pT)d of identified 
7r° mesons as a function of the transverse energy. It has been 
demonstrated earlier that Et is proportional to the total 
number of participating nucleons Npart- For S + Au reactions 
Npart has been calculated for the different centrality cuts by 
Glauber calculations. In Fig.|o|the values of {pT)d as a func- 
tion of Npart are shown. A rise in {pT)d with increasing Npart 
is observed, which again reflects the increase with centrality 
of the slope parameters of the exponential fits in Fig.^ or of 
the local slopes shown in Fig. 0. The values seem to saturate at 
« 290 MeV/c for medium-central to central collisions. These 
relatively small values may be compared with those obtained 
from a complete integration of the fits with Eq. ^ For central 
S + Au collisions one obtains (pt) w 380 MeV/c - these val- 
ues, however, have a large uncertainty due to the low pT ex- 
trapolation. 
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Fig. 10. Truncated mean transverse momentum {pT)d of vr" mesons 
as defined by Eq. ^ plotted as a function of the average number of 
participants Npart- The data points correspond to the 8 Et based 
centrality selections listed in Table|^. The open square shows {pT)d 
for pp collisions extracted from the fit with Eq.Q. A cut parameter 
d — 0.4 GeV/c was used. 



4.2 Comparison to pp and pA data 

The transverse momentum spectra presented for S + Au reac- 
tions unambiguously demonstrate distinct variations as a func- 
tion of event centrality. In order to learn, how these character- 
istic changes can be attributed to differences in the size of the 
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Table 2. Parameters obtained by fitting Eq.gto the invariant cross section of neutral pions from different reaction systems and for different 
classes of centralities. The parameter Ti and T2 are local slopes according to Eq|^ The pp data taken at ^/s = 23.5 GeV have been scaled to 
— 19.4 GeV before fitting according to the recipe given in ||2^. 







Fraction 

of ami 


Po 
(GeV/c) 


n 


Po/n 
(MeV/c) 


Ti at 1 GeV 
(MeV) 


T2 at 3 GeV 
(MeV) 


Xred 


S + Au 


centr. 


7.7% 


5.6 ±0.5 


32.0 ± 2.0 


174.5 ± 2.6 


205.7 ± 1.0 


268.1 ± 3.2 


2.2 


centr. 


25% 


5.7 ±0.3 


32.8 ± 1.2 


174.3 ± 1.5 


204.7 ± 0.5 


265.7 ± 1.9 


3.2 




100% 


6.5 ±0.3 


36.4 ± 1.4 


178.6± 1.4 


206.1 ± 1.8 


261.0± 1.8 


2.2 


per. 


75% 


5.3 ±0.4 


31.6± 2.0 


166.7 ± 2.4 


198.3 ± 0.7 


261.5 ± 3.5 


0.6 


per. 


30% 


4.9 ±0.5 


30.2 ± 2.3 


163.0 ± 2.8 


196.1 ± 0.8 


262.3 ± 4.3 


0.6 


S + S 


centr. 


25% 


5.3 ± 1.4 


31.7± 6.8 


166.2 ± 7.5 


197.7 ± 2.4 


260.8 ± 9.8 


1.0 




100% 


5.6± 1.4 


33.5 ± 6.6 


166.8 ± 6.5 


196.6 ± 8.7 


256.4 ± 8.7 


0.7 


per. 


75% 


6.2 ±4.1 


36.2 ± 20.5 


170.2 ± 11.5 


197.9 ± 2.6 


253.2 ± 19.2 


1.2 


p + p 




100% 


4.9 ±0.3 


34± 1 


146 ± 2 


175 


234 


4.45 



interacting system, it is of interest to compare with data from 
pp and pA collisions in a consistent way. Unfortunately, a rig- 
orous comparison cannot be made free of assumptions, since 
the relevant data from different experiments, in general, do not 
correspond to the same ^/s and rapidity range. The large body 
of experimental data collected for pp reactions, however, may 
be used to extract appropriate scaling laws of transverse mo- 
mentum spectra as a function of ^/s and rapidity y and thus 
allow to adopt these data to the corresponding values of our 
experiment. In this work, we have employed the empirical con- 
cept which is pre sented in Ref . 1 3 1 1 together with data from 
Refs. [|,||j29[||j30|. The minimum bias pp data of Ref. [||^, 
taken at vs = 23.5 GeV and at slightly more forward rapidi- 
ties, were thus scaled appropriately to the conditions of the 
WA80 experiment. The systematic uncertainty introduced by 
this procedure is included in the error bars shown. 

A first comparison of the 7r° spectra from pp collisions with 
that of S-induced reactions was shown in Fig. ^ The slope pa- 
rameters T\oc smoothly grow from approx. 150 MeV at tot ~ 
0.3 GeV/c^ to almost 250 MeV at tot ^ 3.0 GeV/c^, demon- 
strating convincingly the strong concave shape of the trans- 
verse momentum distributions. The local slopes both in pH-p 
and heavy ion data increase in a similar way when going from 
low pt to high pt, thus indicating that the curvature is very 
similar in both casesj^ The main difference is in the systemat- 
ically higher values of the local slopes in heavy ion reactions. 
The slopes of the pT spectra in pp and nucleus-nucleus colli- 
sions differ over the whole px range which indicates the oc- 
curence of strong nuclear effects in particle production. 

The variations seen in the local slopes are reflected in the 
truncated mean transverse momentum (Fig. 10). Here, {pT)d 
has also been calculated for pp data from a parameterization. 
The value of 217 MeV/c is considerably below those for the 



' As discussed above, small differences in the curvature can be seen 
from the fits of Eq. 0. 



S + Au data, even for the most peripheral colUsions. Appar- 
ently, the number of participating nucleons even in the periph- 
eral collisions (w 11) is too large to treat them as identical to 
pp collisions. 

The ratio of multiplicities is shown in Fig.[Tl| as a function 
of Pt for the same set of data. In this presentation, we find 
an increase of the (S + Au)cont/pp ratio with pT by almost a 
factor of 10 over the dynamic range of the experiment and for 
(S + Au)poi/pp an increase by a factor of « 4. We would like 
to point out that this large increase of the cross section ratio at 
a given px may be seen as a consequence of the differences of 
the local slope at lower values of px (see Fig.^. The "high px 
enhancement" of the cross sections therefore appears not to be 
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Fig. 11. Ratio of central (left) and peripheral (right) 200 A-GeY 
S + Au to p + p minimum bias tt" multiplicities as a function of pr- 
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Fig. 12. Transverse momentum dependence of a{pT) from ex- 
pression ^ for the present tt" data from minimum bias S + S and 
S + Au at 200 yl-GeV compared with data from other reaction systems 

an independent feature at high px but may rather be considered 
as the resuh of the different slopes over the complete px range. 

For comparison Fig.|ll| includes estimates of the ratios of 
the number of participants (light grey area at low px) and of the 
number of binary collisions (dark grey area at high px)- The 
pion production seems to scale approximately with the number 
of participants at low pT, but increases much more strongly 
at high pt in central S + Au collisions compared to pp than 
expected from the number of binary collisions. 

A similar variation of the cross section ratios was already 
observed by Cronin et al. in a comparison of p + p and p + A 
data. The systematic behavior could be successfully parameter- 
ized by the phenomenological expression 



E 



(p + A) = A"(^'-)i?^(p + p) 



dp'^ dp^ 



(5) 



where, apart from the expected -^s dependence, the parameter 
oi{pt) was found to depend only weakly on rapidity |35| but 
strongly on pT- At lower px, a turned out to be approximately 
0.7 which in Eq. ^ gives a factor close to A^/^ as expected for 
an opaque target nucleus. For higher transverse momenta, a 
was found to increase with px indicating a growing participa- 
tion of the nuclear volume, approaching a = 1 in the limit of 
full participation of the target nucleus. The surprising obser- 
vation of a > 1 for the highest px values, called 'anomalous 
nuclear enhancement', was subsequently interpreted as an ef- 
fect of multiple scattering of the incident partons . 

For A + B collisions (where A and B denote heavy nuclei), 
the straightforward extension of Eq.|| writes 



B) = {A-B) 



a{pT 



P) (6) 



and has successfully been applied to experimental data in 
Ref. [Q. Fig. 12 shows a compilation of a{px) obtained from 
the present data together with other data sets. The investigated 
projectiles, targets, and detected particle species are indicated 
together with their references. One should note that the statisti- 
cal accuracy of the present data even at high px is comparable 
to the data of [ ^ ] and is superior to that of the other data sets. 



It appears that Eq.^ provides an almost universal parame- 
terization for all systems with a > 1 at high px- The differ- 
ences of the exponent of the present data compared to NA34 
at low px [0 may be caused by the more backward rapidity 
region of the NA34 experiment. However, the data for S H- S 
and S + Au data do also disagree with the data of [ |3^ ] at px — 
0.77GeV/c which are measured in a similar rapidity region. 
It should be noted that the different trigger biases are another 
source of systematic errors in such a comparison of minimum 
bias cross sections. It may also be relevant that the present data 
are compared to a compilation of pp data, while in [B2p one par- 
ticular data set measured within the same experiment has been 
used. Different normalizations of the pp data could therefore 
contribute to a difference in a. 



5 Discussion 

The slope and the curvature of the transverse momentum spec- 
tra, studied as a function of projectile energy, rapidity, and mul- 
tiplicity, carry important information about the reaction dynam- 
ics. All spectra of produced particles emerging from pp, p^, 
and AB collisions are known to exhibit in a reasonable ap- 
proximation an exponential shape which might be regarded as 
indication for a thermal emission mechanism. In fact, for a ther- 
malized system the particle px spectra reflect the temperature 
of the system at the time of decoupling. This idealized picture 
is altered, if there is collective motion in the system, such as 
hydrodynamic flow. In such a case, the transverse velocities of 
particles are to be considered a sum of a random thermal com- 
ponent and an ordered collective expansion velocity. If one or 
both components change between different reaction systems or 
impact parameter regions, the transverse momentum distribu- 
tions are expected to be affected. The picture of thermal emis- 
sion has, despite its simplicity, proven very successful. The as- 
sumptions made in that interpretation might, however, be ques- 
tioned particularly in pp systems because of the small volume 
available for particle production. Surprisingly, the same con- 
cept has recently even been applied to hadron production in 
e+e^ collisions at LEP energies and was able to reproduce 
many of the observed particle ratios [p6[]. One may take the 
success of these models as a hint for the validity of statistical 
concepts, which however do not necessarily imply a thermal 
origin. 

With the growing amount of data and their improved pre- 
cision it has been well established that pT-spectra of all sys- 
tems, when measured with sufficient accuracy and over a large 
Px range, display a curved structure which cannot be fit by a 
simple exponential. Therefore, fits often have been made with 
either a single exponential in a restricted px region and devia- 
tions have then been termed "low px enhancement" or "high 
Px enhancement" respectively, or a superposition of two or 
more exponentials has been applied over a larger range of px- 
Without a physical picture, the former approach is rather mis- 
leading since it arbitrarily defines a "normal" behavior and de- 
viations from that, and the latter one easily leads to an inflation 
of fit parameters. 

On the other hand, there are in fact indications for differ- 
ent mechanisms at high and low px which have to be treated 
separately and should be included appropriately in a complete 



10 



Albrecht et al.: Transverse momentum distributions of neutral pions from nuclear collisions at 200 AGeV 



description of the full px spectrum. For pT ^ 2 GeV/c spectra 
from pp collisions have successfully been described by pertur- 
bative QCD calculations |27|. The "nuclear effects" observed 



in pA reactions have then been discussed in terms of multi- 
ple parton scattering. A generalization and extension of these 
calculations for AB collisions is provided by the kinetic par- 
ton cascade model [ ^7|j3^ ], which treats the reaction - for pro- 
cesses of large momentum transfer - on the basis of interact- 
ing quarks and gluons. Since this model employs perturbative 
QCD, it can only be applied to large px processes. For the 
much more abundant low pT processes, presently only phe- 
nomenological models can be applied. 

Several codes are available, which are based on phenomeno- 
logical string models [|l^Jll|,|l2| and use experimental data from 
elementary collision processes as input. Within these models, 
after appropriately accounting for the nuclear geometry, a heavy 
ion collision is treated as a superposition of individual nucleon- 
nucleon collisions. Obviously, such an approach neglects any 
collective effect which might be present in an extended strongly 
interacting system. Comparison with various data shows that 
such a straightforward extrapolation from pp collisions pro- 
vides a reasonable description, particularly for the intermediate 
Pt range. This indicates that minimum bias heavy ion colli- 
sions in this region do not exhibit dramatic effects extending 
beyond the straightforward extrapolation from pp collisions. 
Deviations from standard string model simulations are observed 
particularly for central nuclear collisions both at low and high 
Pt- Various attempts have thus been made to improve the de- 
scription of the data by including additional mechanisms into 
the string models. These extensions include hard scattering pro- 
cesses [[lip, rescattering of the produced particles in the target 
nucleus and among themselves [p ^ , p roduction of quark clus- 
ters ijs^], or color rope formatioii [[lo||. 

As an example, VENUS 4.12 iQ and FRITIOF 7.02 ^ 
calculations are shown in Fig. [ij together with the 7r° spectrum 
from central Sh-Au collisions. The simulated events have been 
filtered with the experimental acceptance and no normalization 
to the cross section of the data is applied. The description of 
the spectra is reasonably good at Pt < IGcV/c. At higher 
FRITIOF falls below the data and is about a factor of 5 too low 
at > 2 GeV/c. The agreement between data and VENUS 
is better, which may be attributed to the rescattering of pro- 
duced particles implemented in this model. VENUS however 
also starts to deviate at higher pT', it is more than a factor of 
2 too high at Pt > 3 GeV/c. If rescattering in the nuclear 
medium is included, the description at high px is improved 
somewhat, at the expense, however, of additional parameters 
in the model. Still the description is not perfect. 

A completely different approach starts from nuclear hy- 
drodynamics and treats the heavy ion collision as an expand- 
ing thermodynamic al system characterized by a t emper ature, 
chemical potential, and expansion velocity [ pj|jl4l , p4y42| ]. Es- 
pecially for central collisions of heavy nuclei thermal equi- 
librium seems plausible, at least locally, due to the frequent 
rescattering of produced particles in the extended medium. This 
picture is supported by the large rapidity shifts of produced 
particles observed in central compared to peripheral collisions 
1 24 , Jsj] . The extraction of freeze-out temperatures from mea- 
sured pion Pt spectra within this model is complicated for sev- 



eral practical and conceptual reasons. First, it has been recog- 
nized that a significant fraction of the observed pions originate 
from decays of short-lived heavier resonances. This leads to 
characteristic deviations from a single exponential and, partic- 
ularly at low Pt, to a more curved spectrum, as observed ex- 
perimentally [Hp5|. By including all contributing resonances, 
the authors of were able to describe the experimental pT 
spectra of pions from pp systems. The success of this approach, 
however, does not necessarily prove that thermal equilibrium is 
in fact reached in such systems. 

The same procedure has been applied to heavy ion systems 
using a computer program provided by the authors of [p^. 
Here, however, in addition to the most important resonance de- 
cays, transverse flow has also been included in the calculations. 
This leads to a "blue shift" of the apparent temperatures rela- 
tive to the freeze-out temperatures of the system. By fixing the 
baryonic chemical potential to /is = 200 MeV and varying the 
temperature and expansion velocity, good fits to the experimen- 
tal data have been obtained. An example of such a calculation 
applied to the present data is shown in Fig.p^. 

While in earlier analyses of pT spectra these models were 
not sensitive enough to fix simultaneously the values of the 
temperature and the radial expansion velocity without ambigu- 
ities, the precision of the present data allows to further restrict 
the parameter values. The spectrum is best described with a 
temperature T = 184 MeV and an average transverse flow ve- 
locity of (/3) = 0.21. The fits are not very sensitive to the bary- 
onic chemical potential because baryon resonances contribute 
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mostly at very low pT where the experimental coverage van- 
ishes. This can be seen from Fig.n^ where the total yield is 
broken down into the direct thermal pions (dashed line), the 
contribution from p decays (dotted line) and all other reso- 
nances (dash-dotted line). 

Fig. p3| shows 1-sigma-confidence regions in the (T,(/3)) 
parameter plane for these fits. Calculations have been performed 
under the following conditions: 



1 . Including all resonances using 

(a) the full experimental mr-range (solid area in the upper 
part of Fig. |l5|) and 

(b) only data for rriT < 1.6MeV/c^ (outline in the upper 
partof Fig.|l5|) 

and 

2. including only the direct pions using the same rriT regions 
(a) and (W as above (solid area and outline in the lower part 
of Fig. 111). 

A fit to the entire niT range with the full model constrains 
the temperature to the range T = 180 - 190 MeV, a fit to the 
restricted itit range yields essentially the same optimal param- 
eters (T — 178 MeV), but relaxes the constraints - still the 1-a 
range of the temperature is above 160 MeV. The 2-a confidence 
region is essentially unbound towards lower temperatures when 
using only data for ttit < 1-6 McV/c^, one can, however, see 
that such fits for lower temperatures are in contradiction with 
the data at higher ttit- Fig, [iq shows the ratio of fits to the ex- 
perimental data for the optimum value (T = 178 MeV) and 
constrained fits for T = 140 MeV and T = 120 MeV. It can 
be seen that the fits for lower temperatures overpredict the cross 
section at higher uit- In the presence of additional hard scat- 
tering contributions the experimental data would still provide 
an upper bound for the hydrodynamical contribution. There- 
fore the full niT range may be used to estimate a lower bound 
of the temperature parameter. 

It is interesting to note that fits ignoring the resonance con- 
tributions would indicate much lower temperatures and also 
slightly higher flow velocities than those of the full calculation. 
In this case the temperature and flow velocity are rather close 
to the values given in for S H- S collisions (T « 140 MeV, 
{13) = 0.24 - 0.28) Similar numbers ai-e given in [H]. One 
should keep in mind, however, that there is no good reason to 
exclude the resonance decays, and that our results for the full 
calculation are not consistent with the values in [^^Q. 

While the calculation with resonances is clearly more rea- 
sonable than the one without, there are still further uncertain- 
ties related to the model assumptions. Especially the choice of 
one temperature to characterize both chemical freeze-out of all 
resonances and thermal freeze-out of the pions is questionable. 
The chemical freeze-out temperature of the p meson as a res- 
onance appearing in tttt scattering might be closely related to 
the thermal freeze-out temperature of the pions. However, this 
is most likely not true for other particles like e.g. the K*. It 
is therefore necessary to add information from other measure- 
ments, e.g. relative particle abundancies or interferometry mea- 
surements. 

Due to the finite value of the expansion velocity, i.e. an 
energy sharing between random thermal and ordered collec- 
tive motion, the resulting temperature is now smaller than the 
slope parameter obtained from the exponential fits in Fig. || 
The temperature of 184 MeV can be compared to the value of 
T = 160 — 170 MeV extracted from relative particle abundan- 
cies 



The S+S data in the present publication yield similar results to 
the S+Au data, their statistical uncertainty is, however, much larger so 
that they have not been considered further. 
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Fig. 16. Ratio of fits of a thermodynamic model including resonance 
decays Q to data for central (7.7 % of (Jmb) 200 GeV S + Au colli- 
sions for different temperatures. 

6 Summary 

We have discussed transverse momentum spectra of neutral pi- 
ons produced in various systems. The extracted local slopes as 
well as the power-law fits to the pp and AB px spectra un- 
ambiguously demonstrate concave distributions for all systems 
and impact parameter selections. The shapes of the spectra can 
be well characterized by two parameters, e.g. the inverse slope 
at a given px (T = Ti) and the linear variation of the slope pa- 
rameters (AT ^T2- Ti). While AT (the "slope of the local 
slope") and thereby the curvature are very similar in all sys- 
tems, there is a strong increase of the inverse slope parameters 
and of {pt) in going from pp to heavy ion collisions and a con- 
tinuiiig, but more subtle growth with the size of the interaction 
zone. I 

The comparison to event generators shows, that they still 
fail to describe all subtle features observed in the data if mea- 
sured over a large dynamic range. The FRITIOF model shows 
large discrepancies, VENUS 4.12 with its phenomenological 
rescattering provides a better, yet still unsatisfactory descrip- 
tion. A hydrodynamically motivated calculation with two fit- 
ted parameters agrees very nicely with the data. Within such 
a calculation both transverse flow and decay of long-lived res- 
onances appear to be of importance for the interpretation of 
measured transverse momentum spectra. The fit applied to the 
data in the present paper favours comparably large freeze-out 

' Corrections for the influence of a shift of the effective center-of- 
mass on the transverse momentum spectra measured as a function of 
centrality in asymmetric heavy-ion collisions have not been consid- 
ered here, although they might - within certain thermal models - ac- 
count for part of the centrality dependence of the slope parameters. 



temperatures and low flow velocities. As such large tempera- 
tures might require high pion phase space densities, it has to be 
investigated whether they may be interpreted as freeze-out tem- 
peratures or whether the model used here with the assumption 
of local thermal equilibrium may be not applicable. 

The spectra are broadened when going from pp to AB col- 
lisions. At first sight, ratios of the type pA/pp, AB/pp, cen- 
tral/peripheral, etc. seem to isolate the existence of strong high 
Pt phenomena. A closer look at the individual changes of the 
spectral shapes, however, allows to trace the increasing ratios 
with increasing px back to the differences in the overall shape 
of the spectra which necessarily leads to large effects in the 
cross section at large pT- The difference in the slopes over the 
full rriT range may be taken as a hint for non-trivial nuclear ef- 
fects observed in these collisions. The observed variations may, 
in the framework of thermodynamics, be taken as an indica- 
tion of higher temperature and/or a possible stronger collective 
expansion of the system in heavy ion collisions compared to 
pp. Likewise, in the picture of string models, it may suggest a 
stronger amount of rescattering of produced particles. Another 
possible conclusion would be to call for different mechanisms, 
i.e. string formation and decay as sufficient for pp collisions, 
but local thermodynamic equilibrium and hydrodynamic ex- 
pansion in heavy ion collisions. 

Common to all models, even the most conservative ones, 
is the need to either incorporate large amounts of initial- and 
final-state scattering, or to assume high temperature and pres- 
sure build up. Therefore, they imply a dense and strongly inter- 
acting system, i.e. conditions which are favorable to eventually 
reach thermal equilibrium. 
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